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The lightest supersymmetric particle, most likely the lightest neutralino, is one 
of the most prominent particle candidates for cold dark matter (CDM). We show 
that the primordial spectrum of density fluctuations in neutralino CDM has a sharp 
cut-off, induced by two different damping mechanisms. During the kinetic decou¬ 
pling of neutralinos, non-equilibrium processes constitute viscosity effects, which 
damp or even absorb density perturbations in CDM. After the last scattering of 
neutralinos, free streaming induces neutralino flows from overdense to underdense 
regions of space. Both damping mechanisms together define a minimal mass scale 
for perturbations in neutralino CDM, before the inhomogeneities enter the non¬ 
linear epoch of structure formation. We And that the very first gravitationally 
bound neutralino clouds ought to have masses above which is six orders 

of magnitude above the mass of possible axion miniclusters. 


1. Introduction 

Recent measurements support the idea that there is a significant amount of 
cold dark matter (CDM) in the Universe. An analysis of the temperature 
anisotropies in the cosmic microwave background (CMB) gives for the CDM 
mass density flcdinh^ = 0-13io’o2 data with weak /i-prior)^, while from 
the same analysis the mass density of the baryons is much smaller = 
0.023)')q qq 3 . The latter is consistent with the measurement of the primordial 
deuterium abundance in high-redshift hydrogen clouds, which gives the 
most precise determination of the baryon density^, = 0.0205±0.0018. 

Large galaxy redshift surveys support the CMB measurements: from the 
analysis of 160,000 galaxies Percival et al^ find (fib ± f^cdm)^ = 0.20 ± 0.03 
and r2b/(Db ± Dcdm) = 0.15 ± 0.07. Taking various observations together. 
Turner^ estimates that the non-baryonic mass in the Universe, which is just 
the mass of CDM, is given by Dcdm = 0.29 ± 0.04. 
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The characteristic feature of CDM is its non-relativistic equation of 
state at the time when the Universe contains enough matter within one 
Hubble volume to form one galaxy. The two leading particle candidates 
for CDM are the axion, a very light particle with a mass in the range 
Too ~ (10“® — 10“^) eV, which is cold since it has been produced in a 
Bose condensate, and the neutralino, a heavy particle with a mass ~ 
(50 — 500) GeV, which is cold by virtue of its large mass. In the constrained 
minimal supersymmetric extension of the standard model (MSSM), the 
lightest neutralino is most likely the bino. 

For almost 20 years there has been a strong working hypothesis: axion 
CDM and neutralino CDM cannot be distinguished by purely cosmological 
observations because both particle candidates interact only via gravity, as 
far as cosmology goes. This fact has made the study of structure formation 
on large scales > 1 Mpc simple—the microphysics of CDM is irrelevant. 
However, at the galactic scale and below, various particle candidates might 
be distinguishable. To learn more about the nature of CDM, we study the 
small-scale structure of neutralino CDM with emphasis on the very first, 
purely gravitationally bound, neutralino clouds (for details, see our recent 
work®). 


2. Chemical and kinetic decoupling 


There are two distinct temperature scales for CDM particles that are mas¬ 
sive and weakly interacting and obey a thermal history. For temperatures 
T > Ted, binos are kept in chemical equilibrium with all fermions in the 
heat bath via the annihilation processes X + X ^ ^ the rate 
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Here, x = M^/T and Mp denotes the universal sfermion mass. The numer¬ 
ical prefactor is the effective neutralino-fermion coupling, the second factor 
shows the dependence on the MSSM parameters and the x-dependent fac¬ 
tor gives the temperature dependence in units of the neutralino mass and 
is proportional to the bino number density. Chemical decoupling of binos 
happens at the temperature Ted = M^/xed, when the neutralino annihila¬ 
tion rate becomes comparable to the Hubble rate: 
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Scanning the MSSM parameter space, we typically find x^d ~ 22, see Fig¬ 
ure 1. The chemical decoupling temperature Ted increases with increas¬ 
ing bino mass, since the neutralino number density is suppressed with 
exp {—M^/T). As a consequence, the annihilation rate approaches the Hub¬ 
ble rate faster for larger bino masses. For a fixed bino mass, Ted increases for 
an increasing sfermion mass, since the interaction range decreases. (Note 
that this discussion is not correct on a quantitative level, since we neglect 
here important effects like co-annihilation, in order to render an analytic 
discussion possible.) 




Figure 1. Chemical and kinetic decoupling as a function of the universal sfermion (slep- 
ton) mass for bino masses G {50,100,200} GeV (from bottom to top). 


For temperatures Tkd < T < Ted neutralinos are in local thermal equi¬ 
librium (Ite) with all fermions of the heat bath. Equilibrium is maintained 
due to elastic scattering processes x + -P’} ^ X + at the rate 

Fei « 10-2 ^^ x-s . (3) 

Kinetic decoupling of binos takes place when the relaxation time Treiax in 
the bino system exceeds the Hubble time. The relaxation time Treiax = 
Ntcoh differs from the collision time Tcoii = 1/Fei by the number N(T) of 
elastic scatterings needed to keep or establish Ite in the bino system. This 
number is given by the relative momentum transfer per elastic scattering, 
N{T) = {AP^/P^)~^ « M^/T, with AP^ = \/<t> « denoting 
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the rms of the Mandelstam variable t. Kinetic decoupling happens at 
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with a = 0 if mistakenly TcoIi = iVeiax is assumed and a = I when 
the number of scatterings N is taken into account. We typically find 
Tkd = (10— 100) MeV, see Figure 1. The kinetic decoupling temperature is 
increasing with increasing bino mass because the momentum transfer to the 
heat bath is decreasing. As a consequence the number of elastic scattering 
processes needed to keep or establish Ite is increasing. 

We find Ted Tkd, because Fann/Fei « 10~^x®/^ exp (—x) <C 1 for 
T < Ted < the number density of neutralinos is Boltzmann suppressed 
with respect to the number densities of the relativistic fermions. This is the 
reason for the temperature hierarchy Ted Tkd > Tis, with Tis denoting 
the temperature at which binos scatter for the last time. 


3. Local transport coefficients 

During the process of kinetic decoupling, non-equilibrium processes consti¬ 
tute themselves as viscosity phenomena in the bino system. We therefore 
use hydrodynamics for the description. For T ^ Ted, CDM and the heat 
bath can be described by a single ideal fluid. For Ted > T > Tkd the CDM 
fluid is strongly coupled to the radiation fluid (rad), which keeps CDM in 
Ite. Around Tkd, the CDM fluid decouples from the radiation fluid, in which 
Ite persists, since fledm = (u/Geq)llrad ^ Hrad foc T Tgq. 

The resulting non-equilibrium processes in the CDM fluid can be taken 
into account by additional Lorentz tensors and in the current 
density Jedm and the energy momentum tensor Tcdm of the CDM fluid. 
We fix the ambiguities in the relativistic description of imperfect fluids by 
demanding = 0 and = ( hV • U -I- 77 -I- y Sym (U 0 

U denotes the adiabatic velocity field and h projects on the hypersurface 
perpendicular to it. The first term is the bulk viscosity, describing the 
flow of U in the hypersurface defined by h. The second term is the shear 
viscosity and describes the bending of U in the direction perpendicular to 
the adiabatic current. The last term is the heat conduction. 

The strength of the dissipative processes is given by the local transport 
coefficients C, rj and x- An efficient method for calculating these coefficients 
was proposed by Silk®. It has been applied in great detail to the case of a 
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relativistic fluid by Weinberg^. A generalisation to an arbitrary equation 
of state can be found in our recent work®. 

The main idea is to compare the Lorentz tensors in the hydrodynami- 
cal description with the Lorentz tensors in the kinetic description. In the 
kinetic description we make the ansatz {u!,n, x) = A{uj,x) + B (w, x) ■ 
n + C{uj,x) ■ (n(8)n+l/3h) for the CDM phase-space distribution de¬ 
scribing the non-equilibrium state of the bino system. Here, u denotes the 
projection of the momentum on the velocity field and n denotes a vector 
perpendicular to the adiabatic current. We calculate® the coefficients in 
the expansion of into irreducible polynoms in n and h. For the tensor 
structure we find A oc V • V, B oc and C oc The scalar A gen¬ 

erates dissipative processes, which are not perpendicular to the adiabatic 
current. This contribution deserves special care. 

Calculating the energy momentum tensor in the kinetic description 
and comparing it with we find the local transport coefficients C = 

Spcdm/SFei, rj = Pcdm/rei and % = 0 in first order in I/Tei and 1/x. All 
coefficients are decreasing with an increasing elastic scattering rate, since 
elastic scattering allows the transfer of derivations from Ite to the heat bath. 
Heat conduction is a subdominant process for non-relativistic particles and 
vanishes in leading order. 


4. Acoustic absorption 

The dissipative processes presented in the last section transfer energy and 
momentum in the direction perpendicular to the adiabatic flow. This pro¬ 
vides a damping mechanism for acoustic perturbations in CDM^. The 
damping of density inhomogeneities is given by Im(u;) as calculated in lin¬ 
ear perturbation theory for relativistic hydrodynamics. For an acoustic 
wave with wave number k we find 
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where Md « 3 • IO“®(GeV^/M^Tkd)^^^(ffx/i^)Af 0 is the characteristic 
damping scale. We find Md « see Figure 2. Thus, only acoustic 

perturbations with masses M > Md, contained in the overdense volume, 
are not absorbed and enter the free streaming regime. 
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Figure 2. Acoustic damping scale as a function of the bino mass for sfermion masses 
Mp G {150, 200, 300, 400} GeV (from bottom to top). 

5. Free streaming 

For temperatures T < Tig, CDM does freely streaming on geodesics. As 
a consequence, CDM propagates from overdense to underdense regions, 
thus smearing out local inhomogeneities. We find for the induced damping 
scale Mfs(a) = Mdln^(a/ais). Free streaming becomes the leading damping 
mechanism once the universe has doubled its size after last scattering. It 
is interesting to calculate the free streaming scale at the time of matter- 
radiation equality, when CDM density perturbations start to grow linearly 
with the expansion factor. We typically find Mfs(aeq) « 

6. Conclusions 

We have shown that collisional damping and free streaming smear out all 
power of primordial density inhomogeneities in bino CDM below 10 “®Mq 
by the time of matter-radiation equality. This is in striking contrast to 
claims in the literature® that the minimal mass for the first purely gravi¬ 
tationally bound neutralino cloud is 10“^® (150 GeV/M^)® Mq. The huge 
difference to our result stems from the assumption that chemical and ki¬ 
netic decoupling happened simultaneously, which we proved to be wrong. 
We find instead that the very first bino objects have to have masses above 
10“®Mq. This result is very robust with respect to the MSSM parameters. 
These bino clouds are very different from possible axion mini-clusters with 
typical masses® around 

According to hierarchical structure formation these very first CDM 
clouds are supposed to merge and form larger objects. Large scale struc¬ 
ture simulations show structure formation on all accessible scales down to 
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the resolution of the simulation^*^. However, the dynamic range of todays 
simulations is not sufficient to deal with the very first CDM objects, so 
the fate of these objects is an open issue. A cloudy distribution of CDM in 
galactic halos would have important implications for direct^^ and indirect^^ 
searches for dark matter. 
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